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Summary: The Numerical Renormalization Group method (NRG) has 
been developed by Wilson in the 1970's to investigate the Kondo problem. 
The NRG allows the non-perturbative calculation of static and dynamic 
properties for a variety of impurity models. In addition, this method has 
been recently generalized to lattice models within the Dynamical Mean 
Field Theory. This paper gives a brief historical overview of the develop- 
ment of the NRG and discusses its application to the Hubbard model; in 
particular the results for the Mott metal-insulator transition at low tem- 
peratures. 



1 The Numerical Renormalization Group and the Kondo 
problem 

The application of renormalization group (RG) ideas in the physics of condensed 
matter has been strongly influenced by the work of Wilson His 'theory for 
critical phenomena in connection with phase transitions' has been awarded the 
Nobel prize in physics in 1982 This paper deals with one aspect in the 
work of Wilson: the numerical renormalization group (NRG) method for the 
investigation of the Kondo problem. 

The history of the Kondo problem |5| goes back to the 1930's when a resis- 
tance minimum was found at very low temperatures in seemingly pure metals 
Q. This minimum, and the strong increase of the resistance p(T) on further low- 
ering the temperature, has been later found to be caused by magnetic impurities 
(such as iron). Kondo successfully explained the resistance minimum within a 
perturbative calculation for the s-d- (or Kondo-) model ||, a model for magnetic 
impurities in metals. However, Kondo's result implies a divergence of p(T) for 
T — > 0, in contrast to the saturation found experimentally. It became clear that 
this shortcoming is due to the perturbative approach used by Kondo. 

An important step towards a solution of this problem (the 'Kondo problem') 
has been the scaling approach by Anderson B . By successively eliminating high 
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energy states, Anderson showed that the coupling J in the effective low energy 
model diverges. However, the derivation only holds within perturbation theory in 
J and is therefore not necessarily valid in the large J limit. A diverging coupling 
between impurity and conduction electrons corresponds to a perfect screening of 
the impurity spin; the magnetic moment therefore vanishes for T — > and the 
resistivity no longer diverges. This result has been finally verified by Wilson's 
NRG, as will be discussed below. 

In the following, some details of the NRG method are explained in the context 
of the single impurity Anderson model jjj (Wilson originally set up the RG 
transformation for the Kondo model, but the details of the NRG are essentially 
the same for both models |l], |||). The Hamiltonian of this model is given by 

H = ^etftf. + Ufffiflfi 

a 

+5>(# c *-+ c L/*)- (i-i) 



In the model (IT), cj^ denote annihilation (creation) operators for band states 
with spin a and energy Sk, fa those for impurity states with spin a and energy 
£f . The Coulomb interaction for two electrons at the impurity site is given by U 
and both subsystems are coupled via a hybridization V. 
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Figure 1 Logarithmic discretization of the conduction band 
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The first step to set up the RG-transformation is a logarithmic discretization 
of the conduction band (see Fig. [I]): the continuous conduction band is divided 
into (infinitely many) intervals [£ n +i,£n] and [— £ n , — £ n +i] with £ n = DA~ n 
and n = 0,1,2,.... D is the half-bandwidth of the conduction band and A 
the NRG-discretization parameter (typical values used in the calculations are 
A = 1.5, . . . , 2). The conduction band states in each interval are then replaced 
by a single state. Although this approximation by a discrete set of states involves 
some coarse graining at higher energies, it captures arbitralily small energies near 
the Fermi level. 

In a second step, the discrete model is mapped on a semi-infinite chain form 
described by the hamiltonian (see also Fig. 2): 
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Here, the impurity operators are written as f_l a and the conduction band states 

as fnJ with n = 0, 1, 2, . . .. Due to the logarithmic discretization, the hopping 
matrix elements decrease A""/ 2 . This can be easily understood by 

considering a discretized conduction band with a finite number of states M (with 
M even). The lowest energy scale is, according to Fig. [l] given by w Dk~ AI / 2 . 
This discrete model is mapped onto a semi infinite chain with the same number 
of conduction electron degrees of freedom, M. The only way to generate the low 
energy scale ~ DA~ M / 2 is now due to the hopping matrix elements e n so that 
they have to fall of with the square root of Lambda. 
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Figure 2 Semi-infinite chain form of the single impurity Anderson model 

This means that, in going along the chain, the system evolves from high en- 
ergies (given by D and U) to arbitralily low energies (given by DA~ M / 2 ). The 
renormalization group transformation is now set up in the following way. 

We start with the solution of the isolated impurity, that is the knowledge of all 
eigenstates, eigenenergies and matrix elements. The first step of the renormaliza- 
tion group transformation is to add the first conduction electron site, set up the 
hamiltonian matrices for the enhanced Hilbert space, and obtain the information 
for the new eigenstates, eigenenergies and matrix elements by diagonalizing these 
matrices. This procedure is then iterated. An obvious problem occurs after only 
a few steps of the iteration. The Hilbert space grows as A N , which makes it im- 
possible to keep all the states in the calculation. Wilson therefore devised a very 
simple truncation procedure in which only those states (typically a few hundred) 
with the lowest energies are kept. This truncation scheme is very successful but 
relies on the fact that the hopping matrix elements are falling of exponentially. 
High energy states therefore do not change the very low frequency behaviour 
and can be neglected. 
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Figure 3 Flow diagram for the lowest lying energy levels for the single impurity 
Anderson model with e s = -0.2, U = 0.4 and A = 0.015. 



This procedure gives for each cluster a set of eigenenergies and matrix ele- 
ments from which a number of physical properties can be derived (this will be 
illustrated for the calculation of the spectral function in the next section). The 
eigenenergies itself show the essential physics of the Kondo problem: Fig. 3 shows 
the dependence of the lowest lying energy levels on the length of the chain (the 
energies are scaled by a factor A N / 2 ). The system is first approaching an unstable 
fixed point at N m 10 — 20 (the Local Moment fixed point) and is then flowing to 
a stable fixed point for N > 50 (the Strong Coupling fixed point). By analyzing 
the structure of the Strong Coupling fixed point and by calculating perturba- 
tive corrections about it, Wilson (for the Kondo model [uj) and Krishnamurthy, 
Wilkins and Wilson (for the single impurity Anderson model M) found that 

i) right at the fixed point, the impurity spin is completely screened; 

ii) on approaching the fixed point, the thermodynamic properties are Fermi- 
liquid like; i.e. the magnetic susceptibility x{T) approaches a constant 
value for T — ► and the specific heat C — 7T is linear in T for T — ► 0; the 
ratio R — xH ^ s known as the Wilson ratio and takes the universal value 
R = 2 in the Kondo model; 
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2 Developments and Applications of the NRG method 

The NRG approach decribed so far has two main advantages: it is non-pertur- 
bative and can deal with arbitrary values of U (simply because the impurity 
part is diagonalized exactly); and it can describe the physics at arbitrary low 
energies and temperatures (due to the logarithmic discretization). This is impor- 
tant in Wilson's calculation for the Kondo problem which indeed showed what 
had been anticipated by Anderson: the development of a ground state with a 
completely screened impurity (the Fermi- liquid or strong-coupling fixed point). 
The crossover to this fixed point occurs at the Kondo scale 

A\ 1/2 / nU 



k ^ = D {2u) C ^{-Sa)- (2 - 3) 

(This form is valid in the particle-hole symmetric case £/ = —U /2; A is defined as 
A = ^ttV 2 N(Ef) with N(Ep) the density of states of the conduction electrons at 
the Fermi level) . A sufficiently large ratio U / A can therefore generate arbitrarily 
low energy scales. 

On the other hand, the NRG method has one main drawback: it is only ap- 
plicable to impurity type models and therefore lacks the flexibility of e.g. the 
Quantum-Monte-Carlo method. A typical example where the NRG fails is the 
one-dimensional Hubbard model. This model is very similar to the semi-infinite 



chain model of eq. (1.2), but with constant hopping matrix elements between 



neighbouring sites and a Coulomb-repulsion U on each site. One might therefore 



expect a similar iterative diagonalization scheme as for the hamiltonian (1.2) to 
work for the Hubbard model as well. However, the truncation scheme (keeping 
only the lowest lying states) does not work for a model where the same energy 
scales (U and the bandwidth) are added at each step of the RG procedure. The 
low energy spectrum of the cluster with one additional site now depends on states 
from the whole spectrum of energies of the previous iteration. (A solution to this 
problem, i.e. finding a truncation scheme which gives an accurate description of 
the larger cluster, is the Density matrix renormalization group method |^]). 

There are, fortunately, a lot of interesting impurity models where the NRG 
can be applied and where it provided insights into a variety of physical problems. 
Non- Fermi liquid behaviour has been studied in the context of the Two-Channel- 
Kondo-Model and related models |l(J. The structure of the Non- Fermi liquid 
fixed point as well as its stability against various perturbations has been clarified 
using the NRG method. 

Another example is the quantum phase transition in impurity models coupling 
to conduction electrons with a vanishing density of states at the Fermi level: 
p c (w) <x |w| r . Here the NRG enables a non-perturbative investigation of both 
the strong-coupling and local moment phases as well as the quantum critical 
point seperating these two JTTI] . 



G 



Bulla 



Apart from applying the NRG to generalized impurity models, some impor- 
tant technical developments have been made during the past 10 - 15 years; most 
notably the calculation of dynamical properties, both at zero and finite temper- 
atures n || . 

Let us briefly discuss how to calculate the single-particle spectral function 
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A{u) = hnG(u + i6+) , with G{z) = ((/,,/])), 



(2.4) 



within the NRG approach. Due to the discreteness of the Hamiltonian, the spec- 
tral function A{u>) is given by a discrete set of <5-peaks and the general expression 
for finite temperature reads: 



A N {w) 



-E 

Z N Z - 



J-lc 



- (E n - E m )) 



-pE„ 



+ e-^«) . (2.5) 



The index N specifies the iteration number (the cluster size) and for each N the 
spectral function is calculated from the matrix elements (n|/i lo .|rn^ and the 



eigenenergies E n ,E m . Zn is the grand canonical partition function. Eq. (2.5) 



defines the spectral function for each cluster and a typical result is shown in Fig. 



Here, the weight of the <5-peaks in eq. (2J5) is represented by the height of the 
spikes. One can clearly see the typical three peak structure from the result of the 
14-site cluster: charge fluctuation peaks centered atww ±0.7 (w ±U /2) and a 
quasiparticle peak at the Fermi level (here u = 0). However, the resolution of the 
quasiparticle peak appears to be rather unsatisfactory: there is no information 
on the spectral density below |w| « 0.04. The advantage of the NRG approach 
(as compared to e.g. the Exact Diagonalization technique) is that by successively 
increasing the length of the chain, one can extract the information on the spectral 
density down to arbitrarily low energy scales. This is seen in the results for the 
N = 16 and N — 18 clusters in Fig. [|. The necessary truncation of states, 
as decribed in the previous section, is also obvious from Fig. ^. There are no 
excitations for |u>| > 0.85 (|w| > 0.45) in the N = 16 (N = 18) cluster, so that 
the information on the charge fluctuation peaks is lost for the N = 16 and larger 
clusters. In order to obtain the spectral density for all energy scales, the data 
from all cluster sizes have to be put together. This means that each cluster size 
only provides the information on its relevant energy scale. 

The resulting spectrum will still be discrete, of course, with the 5-peaks getting 
closer and closer together for to — > 0. It is convenient (both for using the results 
in further calculations and for visualizing the distribution of spectral weight) to 
broaden the 5-peaks in eq. (p.5|) via 



6{lu 



b0Jr. 



exp 



(lnw — lncj„) z 
62 



(2.6) 



The Numerical Renormalization Group Method for correlated electrons 



7 



N=14 



0.20 



0.10 



o.oo - — ^ — LJ.l L.I i. i i J . J ,j, J — Li_. . 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

(0 

N=16 



0.20 



0.10 



0.00 



jJiLl 



llilljj 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

to 



N=18 



0.20 



0.10 



0.00 




-1.5 



-1.0 



-0.5 0.0 



0.5 



1.0 



1.5 



Figure 4 Spectral functions Ajv(w) for clusters with size N = 14, 16 and 18. 
The weight of the 5-peaks is given by the height of the spikes. 
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Figure 5 Spectral functions for the single impurity Anderson model for various 
values of U/A. 



The broadening function is a gaussian on a logarithmic scale with width b. In 
this way, the broadening takes into account the logarithmic distribution of the 
(5-peaks. 

Typical results for the spectral function of the single impurity Anderson model 
are shown in Fig. |^. The spectra clearly show the narrowing of the quasiparti- 
cle resonance on increasing the ratio U/A - corresponding to the exponential 
dependence of the low energy scale 7k on U/A. 

Let us now discuss another, very important development which made it possi- 
ble to apply the NRG method also to lattice models of correlated electrons: the 
Dynamical Mean Field Theory (DMFT). 

Metzner and Vollhardt ju) showed that one can define a non-trivial limit 
of infinite spatial dimensions for lattice fermion models (such as the Hubbard 
model). In this limit, the self energy becomes purely local which allows the 
mapping of the lattice model onto an effective single imp urity Anderson model. 
This impurity model has the same structure as in eq. ( |l.l| ), but the density 
of states of the conduction band in the impurity Anderson model has to be 
determined self-consistently and therefore acquires some frequency dependence. 
The NRG can nevertheless be applied to this case (for details see ) . The first 
attempts to study the Hubbard model is the work of Sakai and Kuramoto flqj . 
The results obtained later by Bulla, Hewson and Pruschke [H and Bulla flTfl 
will be discussed in the following section. 
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3 NRG results for the Mott-Hubbard metal-insulator tran- 
sition 

The Mott-Hubbard metal-insulator transition []l8|, is one of the most fasci- 
nating phenomena of strongly correlated electron systems. This transition from 
a paramagnetic metal to a paramagnetic insulator is found in various transi- 
tion metal oxides, such as V2O3 doped with Cr [ pp| . The mechanism driving 
the Mott-Hubbard transition is believed to be the local Coulomb repulsion U 
between electrons on a same lattice site, although the details of the transition 
should also be influenced by lattice degrees of freedom. Therefore, the simplest 
model to investigate the correlation driven metal-insulator transition is the Hub- 
bard model 

H = -t ( c L c j<?+ c ]° Ci ^ + U J2 c h c ^4i c il> ( 3 - 7 ) 

<ij>a i 

where cj_ (ci a ) denote creation (annihilation) operators for a fermion on site 
i, t is the hopping matrix element and the sum Yl<ij> is restricted to nearest 
neighbors. Despite its simple structure, the solution of this model turns out 
to be an extremely difficult many-body problem. The situation is particularly 
complicated near the metal-insulator transition where U and the bandwidth 
are roughly of the same order and perturbative schemes (in U or i) are not 
applicable. 

The DMFT has already been briefly discribed in section 2; this method enabled 
a very detailed analysis of the phase diagram of the infinite-dimensional Hubbard 
model p4],|2jj. The nature of the Mott-transition, however, has been the subject 
of a lively debate over the past five years (see ||(], |7[ ||FJ). This debate 

focusses on the existence (or non-existence) of a hysteresis region at very low 
temperatures. In such a region, two stable solutions of the DMFT equations 
should exist: a metallic and an insulating one. This scenario has been proposed 
by Georges et al. based on calculations using the Iterated Perturbation Theory 
(IPT), Quantum Monte Carlo and Exact Diagonalization p4[ . The validity of 
this result has been questioned by various authors p7j , p8| , |29| . 

Let us now discuss the NRG results for the infinite dimcnional Hubbard model, 
first of all for T = 0. The spectral function A(uS) for the Bcthc lattice is shown 
in Fig. I for U = 0.8U C , U = 0.99U C and U=1.1U C (U c » 1.47W, W: bandwidth) 
In the metallic phase (for large enough values of U) the spectral function shows 
the typical three-peak structure with upper and lower Hubbard bands centered 
at ±[7/2 and a quasiparticle peak at the Fermi level. For U = 0.99U C , the quasi- 
particle peak in both Bethe and hypercubic lattice seems to be isolated (within 
the numerical accuracy) from the upper and lower Hubbard bands, similar to 
what has been observed in the IPT calculations for the Bethe lattice . Conse- 
quently, the gap appears to open discontinuously at the critical U (whether the 
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Figure 6 Spectral functions for Bethe lattice for various values of U. A narrow 
quasiparticle peak develops at the Fermi level which vanishes at the critical [/ c ~ 
1A7W. 

spectral weight between the Hubbard bands and the quasiparticle peak is exactly 
zero or very small but finite cannot be decided with the numerical approach used 
here). 

The quasiparticle peak vanishes at U c « 1A7W in excellent agreement with 
the results from the Projective Self-consistent Method (PSCM) §|, |lj U c sa 
1.46M^. Coexistence of metallic and insulating solutions in an interval U c ,i < 
U < U c ,2 is also found within the NRG approach. Starting from {7 = 0, the 
metal to insulator transition occurs at the critical C/ c ,2 with the vanishing of 
the quasiparticle peak. Starting from the insulating side, the insulator to metal 
transition happens at U Ct \ < U c< 2 (the NRG and IPT give U c ,i ~ 1.25W for the 
Bethe lattice). 

The NRG method for the Hubbard model has only recently been generalized at 
finite temperatures |32| . Preliminary results for the spectral function are shown 
in Fig. |^ for T = 0.00625VF and increasing values of U. The upper critical U is 
given by U Cy 2 ~ 1.24W and the transition at {J c 2 is of first order, i.e. associated 
with a transfer of spectral weight. The 'insulator' for U > U c ^ does not develop 
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Figure 7 Spectral functions for T = 0.00625M 7 and increasing values of U. 



a full gap (this is only possible for U — > oo or T — > 0), but the corresponding 
transport properties in this temperature range will be certainly insulating-like. 

For this temperature, the NRG again finds two stable solutions in an inter- 
val U c ,i(T) < U < U C> 2(T); a metallic one, with a quasiparticle peak at the 
Fermi level and an 'insulating' one, with very small spectral weight at the Fermi 
level (not shown here). The exact shape of the hysteresis region has still to be 
determined and will be discussed elsewhere Q . 

What we have seen in this section is that the NRG-method (together with the 
DMFT) can be applied to the infinite-dimensional Hubbard model and allows 
a non-perturbative calculation of dynamical properties. The calculations can be 
performed for arbitrary interaction strength and temperature, so that the phase 
diagram can be (in principle) determined in the full parameter space. 



4 Further developments of the NRG method 

As we have discussed in the previous sections, Wilson's NRG can be applied to 
two different classes of problems: impurity models and lattice models (the latter 
ones, however, only within the DMFT). 

Concerning impurity models, the NRG has provided important theoretical 
insight for a variety of problems and certainly will do so in the future. In the light 
of the increasing possibilities of experimental fabrication, new classes of impurity 
models are becoming of interest. The behaviour of electrons in quantum dots, 
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for example, can be interpreted as that of an impurity in a conduction band 
(for an application of the NRG method to this problem, see jj3)). Magnetic 
impurities can also serve as sensors, put into certain materials in a controlled 
way. Here one might think of impurities in a correlated host [Q, or impurities 
in a superconducting or magnetic medium. A lot of theoretical work in applying 
the NRG method to these problems still needs to be done. 

The second class of models are lattice models within the DMFT. Here, the 
NRG allows (at least in principle) the calculation of a large set of experimentally 
relevant quantities for a wide range of parameters (especially low temperatures 
and strong correlation) for a large class of models. Apart from the application to 
the Hubbard model which has been briefly discussed in section 3, the NRG has 
already been applied to the periodic Anderson model |H| and to the problem 
of charge ordering in the extended Hubbard model |56| . Future work will focus 
on generalizing the NRG method to magnetically ordered states and to systems 
with a coupling to (dynamical) phonons. 

Of particular interest is the generalization of the NRG to multi-band models. 
In this way, the NRG could further extend the range of applicability of the 
LDA+DMFT approach |3?|] . Here, the non-interacting electronic band structure 
as calculated by the local density approximation is taken as a starting point, 
with the missing correlations introduces via the DMFT. On a more fundamental 
level, the basic physics of multi-band models at low temperatures still needs to 
be clarified, and again, the NRG is the obvious choice for investigating such 
models in the low T and intermediate to large U regime. 

The author would like to thank T. Costi, D.E. Logan, A.C. Hewson, W. Hof- 
stetter, M. Potthoff, Th. Pruschke, and D. Vollhardt for stimulating discussions 
and collaboration over the past few years. Part of this work was supported by 
the Deutsche Forschungsgemeinschaft, grant No. Bu965-1/1 and by the Sonder- 
forschungsbereich 484. 
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